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Abstract: We consider a modification of the three-dimensional Navier-Stokes equations and other hydrodynamical 
evolution equations with space-periodic initial conditions in which the usual Laplacian of the dissipation operator is 
replaced by an operator whose Fourier symbol grows exponentially as e^'^'^*'' at high wavenumbers \k\. Using estimates 
in suitable classes of analytic functions, we show that the solutions with initially finite energy become immediately 
entire in the space variables and that the Fourier coefficients decay faster than e"'-'*^^*''^'"'''^!^*''^ for any C < 1/(2 In 2). 
The same result holds for the one-dimensional Burgers equation with exponential dissipation but can be improved: 
heuristic arguments and very precise simulations, analyzed by the method of asymptotic extrapolation of van der 
Hoeven, indicate that the leading-order asymptotics is precisely of the above form with C = C* = 1 / In 2. The same 
behavior with a universal constant C* is conjectured for the Navier-Stokes equations with exponential dissipation in 
any space dimension. This universality prevents the strong growth of intermittency in the far dissipation range which 
is obtained for ordinary Navier-Stokes turbulence. Possible applications to improved spectral simulations are briefly 
discussed. 



1. Introduction 

More than a quarter of a millenium after the introduction by Leonhard Euler of the equations of incompressible fluid 
dynamics the question of their well-posedness in three dimensions (3D) with sufficiently smooth initial data is still 
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moot |[Tl|2][3]|4l (see also many papers in (5] and references therein). Even more vexing is the fact that switching to 
viscous flow for the solution of the Navier-Stokes equations (NSE) barely improves the situation in 3D ll6H7l[8ll9l fT0ll . 
Finite-time blow up of the solution to the NSE can thus not be ruled out, but there is no numerical evidence that this 
happens. 

In contrast, there is strong numerical evidence that for analytic spatially periodic initial data both the 3D Euler and 
NSE have complex space singularities. Indeed, when such equations are solved by (pseudo-)spectral techniques the 
Fourier transforms of the solution display an exponential decrease at high wavenumbers, which is a signature of com- 
plex singularities fTTII . This behavior was akeady conjectured by von Neumann ||T21 who pointed out on p. 46 1 that the 
solution should be analytic with an exponentially decreasing spectrum. Recently Li and Sinai used a Renormalization 
Group method to prove that for certain complex-valued initial data the 3D NSE display finite-time blow up in the real 
domain (and, as a trivial corollary, also in the complex domain) lfT3l . 

For some PDEs in lower space dimensions explicit information about the position and type of complex singularities 
may be available. For example, complex singularities can sometimes be related to poles of elliptic functions in connec- 
tion with the reaction diffusion equation [ 14] and 2D incompressible Euler equations in Lagrangian coordinates ifTSl . 
The best understood case is that of the ID Burgers equation with ordinary (Laplacian) dissipation^ its singularities 
are poles located at the zeroes of the solutions of the heat equation to which it can be mapped by using the Hopf-Cole 
transformation (see, e.g., I116II17I and references therein). 

We now return to the 3D NSE with real analytic data. It is known that blow up in the real domain can be avoided 
altogether by modifying the dissipative operator, whose Fourier-space symbol is yu|^P, to a higher power of the Lapla- 
cian with symbol i-i\k\^" {a > 5/4) ||6l [T8l . The numerical evidence is however that complex singularities cannot be 
avoided by this "hyperviscous" procedure, frequently used in geophysical simulations (see, for example, |fT9l ). 

Actually, we are unaware of any instance of a nonlinear space-time PDF, with the property that the Cauchy problem 
is well posed in the complex space domain for at least some time and which is guaranteed never to have any complex- 
space singularities at a finite distance from the real domain. In other words the solution stays or becomes entire for all 
f > 0. Here we shall show that solutions of the Cauchy problem are entire for a fairly large class of pseudo-differential 

' The case of the Burgers equation with modified dissipation will be considered in Section|5] 
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nonlinear equations, encompassing variants of the 3D NSE, which possess "exponential dissipation", that is dissipation 
with a symbol growing exponentially as e'*'^'^'' with the ratio of the wavenumber \k\ to a reference wavenumber A;d. 

The paper is organized as follows. In Section |2] we consider the forced 3D incompressible NSE in a periodic 
domain with exponential dissipation. The initial conditions are assumed just to have finite energy. The main theorem 
is established using classes of analytic functions whose norms contain exponentially growing weights in the Fourier 
space II20II2TI . In Section [3] we show that the Fourier transform of the solution decays at high wavenumbers faster 
than exp(-C^lnA^) for any C < l/(21n2). Here, k := \k\/kd is the nondimensionalised wavenumber. In Section|4]we 
briefly present extensions of the result to other instances: different space dimensions and dissipation rates, problems 
formulated in the whole space and on a sphere, and different equations. In Section|5]we then turn to the ID Burgers 
equation with a dissipation growing exponentially at high wavenumbers, for which the same bounds hold as for the 
3D Navier-Stokes case. However in the Burgers case, simple heuristic considerations (Section lSTt and very accurate 
numerical simulations performed by two different techniques (Sections 15 .21 and | 5 . 3l l. indicate that the leading-order 
asymptotic decay is precisely exp ((- 1 / In 2)k In k). We observe that the heuristic approach, which involves a dominant 
balance argument applied in spatial Fourier space, is also applicable to the 3D Navier-Stokes case with exactly the 
same prediction regarding the asymptotic decay. In the concluding Section|6]we discuss open problems and a possible 
application. 

2. Proof that the solution is entire 



We consider the following 3D spatially periodic Navier- 

du 

dt ^ 
u(x, 0) = uq(x). 



-Stokes equations with an exponential dissipation (expNSE) 

-fiDu + f, V-M = 0, (2.1) 

(2.2) 



Here, D is the (pseudo-differential) operator whose Fourier space symbol is e ' that is a dissipation rate varying 
exponentially with the wavenumber \k\, uo is the initial condition, / is a prescribed driving force and fi and cr are 
prescribed positive coefficients. The problem is formulated in a periodic domain Q (for simplicity of notation we take 
Q - [0, 2n]^). The driving force is assumed to be a divergence-free trigonometric polynomial in the spatial coordinates. 
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For technical convenience we use cr in the statements and proofs of mathematical results, while the use of the reference 
wavenumber = l/(2cr) is preferred when discussing the results. 

The initial condition is taken to be a divergence-free periodic vector field with a finite norm (finite energy). 

As usual the problem is rewritten as an abstract ordinary differential equation in a suitable function space, namely 

— + H&^'^^"~u + B{u,u)^ f; (2.3) 
dt 

m(0) = Mo; (2.4) 

where A :- -V^ and B{u, u) is a suitable quadratic form which takes into account the nonlinear term, the pressure term 
and the incompressibility constraint (see, e.g. |l6]|2l|9l). Note that the Fourier symbol of A'^^ is |^|. 

The problem is formulated in the space H :- {(f e (L^(Q))^ : ip is periodic, J ipdx - 0,^ ■ ip - 0). Here, for any 
A>0, the Fourier symbol of the operator e''^"'is given by e'*'*', where k e Z^\{(0, 0, 0)). 

To prove the entire character, with respect to the spatial variables, of the solution u{t) of expNSE for f > 0, it suffices 
to show that its Fourier coefficients decrease faster than exponentially with the wavenumber |^|. This will be done by 
showing that, for any A > 0, the L? norm of e'^"'u, the solution with an exponential weight in Fourier space, is finite. 
As usual, we here denote the L? norm of a real space-periodic function / by |/| := ,^^^^^^^^7i/WP^. Moreover, H"' 
will be the usual Sobolev space of index m (i.e., functions which have up to m space derivatives in L^). 



The main result (Theorem 2.1) will make use of the following Proposition which was inspired by 1201 (see also 
Proposition 2.1 Let a > , /3 > 0, k > 3 and ip e dom (e(«+«'*"'). Then 



y'"'B(^,j <c^(mf^''k'"''pr''' |e(-/^)^">^*^ (2.5) 



where Ca is a universal constant and 



1,{J3) sup x^Q-l^-' = 1^1 e"''. (2.6) 



0<:t<oo 
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Notation In Proposition 2. 1 and also in the sequel we use the following notation (to avoid fractions in exponents): 

5 2k + 5 

e„»li_|, ,2.7) 

Proof Let w e H. By using the Fourier representations (f{x) - Yji^^'^'fi ^"d w{x) - Ji^: e'^ '^w^ and Parseval's 
theorem, we have 



1 ( 



(2.8) 



where the * means complex conjugation. 

Since e"'*' < when k = I + m, we can estimate the absolute value of the right-hand side from above as 



<2 Z 

ki^Q l+m=k. l,mi=0 

where the functions 0(x), 'Fix) and W(x) are given by 



(2.9) 



0(x)¥(x)W(x)dx < ■j:^\mL-' 1^1 \W\ , 



1 

(2ny 



and 



0(x)^Yj^"%,\e'-\ (2.10) 

?'(;.) = 2 |m|e«l'"l|^„|e"" ^ (2.11) 



W(x) = Yj^Wk\s''\ (2.12) 



ki^Q 

and the last inequality follows from the Cauchy-Schwarz inequality. 
By Agmon's inequality 1221 (see also Q) in 3D we have 

II'^'IIl" < CaII'/'II^JI'/'II^, < Ca\A""P\^ \A¥\'2 = Ca Iac-^^''^! Ia^^V-^''^! , (2.13) 

where Ca > is a universal constant. By using (12.101 1. (12.131 1 and the fact that |W| = we obtain 

|(e"^"'B((/j,(^),w)| < CAle^^'^J Ue'"^"V|" Uh^^'^J" \w\. (2.14) 



6 Claude Bardos, Uriel Frisch, Walter Pauls, i 'li mriddhi Sankar Ray, and Edriss S. Titi 

And by using the interpolation inequality between L? and H", where a- > 3, we obtaiq^ 

|(e"^^B(^o,^),w)|<CA|e'"*^^>| pe'"^^o| |w|. (2.15) 
Now, to obtain the inequality in Proposition 2.1, we just need to estimate the X.(H) operator norm 

|Lfe-M'/^|| < sup x'^e-'^*' = (-) e-'' = Z,06). (2.16) 

" ''-C{H) 0<x<co 

This concludes the proof of Proposition 2. 1 . 

Next, we state and present the proof of the main theorem. The steps of the proof are made in a formal way, how- 
ever, they can be justified rigorously by establishing them first for a Galerkin approximation system and using the 
usual Aubin compactness theorem to pass to the limit (see, e.g. ll6]|7]|9]). Furthermore, we do not assume that the 
initial condition mq is entire; it is only assumed to be square integrable, although it will become entire for any f > 0. 
This is why in estimating norms of the solution with exponential weights we have to stay clear of f = 0. 

Theorem 2.1 Let uq e H, fix T > and let f{t) — /(. , f) be an entire function with respect to the spatial variable 
X. Then for every n = 0, 1 , 2 . . . there exist constants C„,C„, K„ and K„ which depend on |mo|, yU, T, cr and on the norm 



f 

Jo 



T . 2 

e<"-"'"'^"7(i) ds, (2.17) 



moreover there exists integers Pn,1n^ 1 such that 



1 K 

|e'"^'^"'M(f)| <-^+C„, forall te(0,T] (2.18) 



and 



f \e^"^^'>''^"\(s)\^ds < ^ + C„, forallt e(0,T], (2.19) 

where u(t) is the solution of ( |2.3l )-( |2!4l l. 

Corollary 2.1 Let uo e H,T > and let f{x, t) be an entire function with respect to the spatial variable x such 

T I 1/^ \^ 

that for every M > we have e'^'"'*"V(i) ds < oo. Then, the solution u(t) of ( I2.3b -( I2.4| | is an entire function with 
respect to the spatial variable for all t e (0, T], and satisfies the estimates (|2.18l l and ( |2.19l l in Theorem 2.1 for any 



The simplest formulation is obtained for k = 5 but the optimization of the bound for the law of decay in Section[3]requires using arbitrary k. 
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Proof of Corollary 2. 1 Consider the Fourier series representation 

M(x,f) = ^e''^-'M(;t,f). (2.20) 

k 

From (12.181 ) and Parseval's theorem, we have, for any n = 0, 1, 2 . . . 

2e2'"^l'^l|M(/t,f)P < oo, (2.21) 

k 

for f > 0. In ( |2.20t we change x to a complex location z-x + \y and obtain 

u{x + iy, f) = ^ Q^H'^+'^y) f) (2.22) 

k 

^ J^[e'^.^ e-W] [e-^ >'^l*l M(/t, f)] ■ (2.23) 

k 

For any n- 1,2,..., the series (|2.23l l of complex analytic functions converges uniformly in the strip |y| + 1 < crn. This 
is because the sum in (12.231 1 is shown to be bounded, for any y, by use of the Cauchy-Schwarz inequality applied to 
the two bracketed expressions and use of ( |2.21| i with \y\ + \ < crn. Hence the Fourier series representation converges in 
the whole complex domain. This concludes the proof of the entire character of the solution with respect to the spatial 
variables. 

Remark This corollary just expresses the most obvious part of the Paley- Wiener Theorem. 



P roof of Theorem 2.1 The proof of the theorem proceeds by mathematical induction. 

Step n = We prove the statement of the theorem for n = 0. We take the inner product of ( |2.3l l with u and use the 
fact that (Z?(m, m), m) = to obtain (when there is no ambiguity we shall henceforth frequently denote u{t) by u) 

1 ^|„|2 + |e-A"^/ ^ ^ (e-^"7,e'^'*"'«) (2.24) 

2 at I I 

< e / e'^'* M 

<^-^ + f|e-"^«|^ (2.25) 
2n 2 

where Young's inequality has been used to obtain the third line. Therefore 

|e-^"'«|' < iL!^. (2.26) 
at \ \ u 
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Integrating the above from to T, we obtain 



- ;u f \&''^'^\{sf ds < Co := |moP + - f 
Jo ' I Jo 



M(f)r + e"^-^ u(s)\ ds < Co := M + - k''^ /W ds. (2.27) 



// f |e'^^"'M(i)| ds<Co. (2.29) 
Jo ' ' 



Hence 

\u(tt < Co, (2.28) 

and 

JO 

From (12.281 1 and (12.29b we obtain ( 12.18b and ( 12.19b for the case n = 0. Here Co is given by ( 12.27b . /To = 0, A'o = and 
Co - Co- Notice that since Ko - Ko = there is no need to determine the integers po and ; however, for the sake of 
initializing the induction process we chose po - qo - 1 • 

Step n n + \ Assume that ( 12.181 ) and ( 12.191 ) are true up to n = m and we would like to prove them for n - m + 1. 
Let us take the inner product of (12.31 ) with e^^'"^'^°"'*"^M and obtain 

2dt\ 1^1 I 

< |(/,e2('«+')-'^"'M)| + |(B(M,M),e2('"+i)-^"'M)| 

<|e'"'^^"V| |e('«+2K^"'M| + |e'"'"^"'B(M,M)| . 
Now we use Proposition 2. 1 to majorize the previous expression by 

< je^-^'Vl + Ca (/.(;e))^« |e"-^"'«|'"'*'^ |e("™"'Mp*'^ . (2.30) 

By Young's inequality we have 

CA[/.(yS)]^« le-'-^'^^f |e('-2>-^'%p*^' < 



2 



(2.31) 



It follows that 



dt ' II I ju I I 

+ ?^^-dWc2'W (/,(^))2-WfW (^^f'' U.A>^|2«W 
2/<: \ K \ 



(2.32) 
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Now we integrate this inequality on the interval {s, t) c (0, T), obtaining 

|e<'"+i>-^"'M(f)f +;« £ |e<'«+2)-^"'«(i')f ds' 

< \(m+\)<TA'l^ 



u{s)f + - {' f{s')f ds' + C r'|e'"'"^"'M(s')f'^'' t// (2.33) 

< |e("-"-^"'«(.)f + - C |e''-^"V(i')r ds' + C f |e"-^"'«(*')r'*'^ ds' , 
I I ^ Jo ' ' ' ' 



where we have set for brevity 



C = C(;., A K) := ^ ;u-^« Cf (/,Ce))2^«'M j , (2.34) 

and where l^ip) is given by ( I2.6I I. 

Now we come to the point where we use the actual induction assumptions. We use ( 12.181) and the midpoint convexity 
to estimate the integrand in the last integral: 

|e'"-^"'«(f)| < + C,„) <2'«(^) +2'WC^W. (2.35) 

Whence it follows that 

1 



C f |e"-^"'«(/)|''*'^ ds' < I'^'^C — i ( —\ 

Js 1 1 e(^)p,„-l 



Discarding the positive term fi e^'"^^^'^'^^ '^(^0 ds' in (I2.33K we obtain from < |2.33t and ( |2.361 t 

|e('"+i'-^'%(0|' < |e<'"+i>-^"'«(^)|' + - r |e'«-^"V(/)P ^/i' 



(2.36) 



Jo 

• 2'<'^C' K^M—L^ + 2'Wc' cy(t - s). 

e(K)p,„ - 1 "' 5eW;)„,-I 



(2.37) 



Integrating this inequality with respect to s over (f/2, f) we get 

^(m+l)o-A"- 



e' 



«(f)r ^ - r ie<'""'>'^^"'«(*)r + - r ie'"'^^"V(*')r 

' t J, 12 I ' Jo ' ' 



(2.38) 



(e(/(:)/?,„ - l)(e(/f)/?,„ - 2) \ f ^ 
Note that /9,„ > 1 implies that 

e{K)p,n - 2 > 0. (2.39) 
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By using ( |2.19l l. we have 



2 2 2 r I „, „|2 



|e^'"+"-^"-M(f)| <K,„\-\ +C„,- + - f |e"--^"7(^')| 



(eW/7„-l)(eWp„-2)\^/ " 4 

From this relation follows that ( 12.181 1 holds for m + 1 with 



By the induction assumption we use ( 12.33b to estimate 



J//2 



„(m+I)o-A'' 



(2.40) 



/7,„+i = max{e(/f)/?,„ - 1, + 1) , (2.41) 
q,n+i = max{e(Af)/?„, - 1, + 1) . (2.42) 



• Cm. (2.43) 



From this estimate and the above we conclude the existence of the constants K,„+\,Cm+\ and the integer /?„,+i such that 
( |2.18t holds for n - m + I. Using the estimate that we have just established in ( 12.181 ) for n - m + \, and substituting 
this in ( I2.33I 1. we immediately obtain the estimate ( 12.191 ) for n - m + \. This concludes the proof of Theorem 2.1. 



3. Rate of decay of the Fourier coefficients 

The purpose of this section is to specify the behavior of various constants appearing in the preceding section to 
obtain the rate of decay with the wavenumber of the Fourier coefficients u(k, t) for f > 0. We again consider the 
3D case in the periodic domain. Since the decay may depend on the rate of decay of the Fourier transform of the 
forcing term f{x, f), for simplicity we assume zero external forcing, which we expect to behave as the case with 
sufficiently rapidly decaying forcing. The adaptation to sufficiently regular forced cases, for example a trigonometric 
polynomial, is similar but more technicaiy Furthermore, it is enough to prove the decay result up to a time T such that 
1/A^ TU/L < 1, where L and f/ are a typical length scale and velocity of the initial data. Extending the results to 
later times is easy (by propagation of regularity). 

It is conceivable that the results can be extended to forces entire in the space variables whose Fourier transforms decrease faster than e"'''*' 
with sufficiently large C. 
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We shall show that the bound for the square of the L? norm of the velocity weighted by t""'^ is a double expo- 
nential in n. Specifically, we have 



Theorem 3.1 Let u{t) be the solution of ( I2.3l l-( I2.4I ) in [0, T] with f - and < T < L/U. Then for every k > 3 
and 6 > 0, there exists a number A, depending on 6 and k, such that, for all integer n > 

c"^^"'u(t)f <(^]"" , te{0,T], (3.1) 



|e(«+'>-^"'M(.)|' ds < (^)"" , f e (0, T] (3.2) 

where a„ = ((1 +6) -\ ■ (3.3) 

\ 2k -5 



Corollary 3.1 For any t > the function u(t) of ( 12. 3t - ( 12741 ) is an entire function in the space variable and its 
( spatial) Fourier coefficients tend to zero in the following faster-than-exponential way: there exists a constant A such 
that, for any Q < s < \, we have 



\u(k,t)\ < e for all \k\ > 



[al 




Ut 





(3.4) 



where 



4-K-5 
'2/C-5 



(3.5) 



Proof of Corollary 3.1 Since we are dealing with a Fourier series, the modulus of any Fourier coefficient of the 
function e*"^''°"'* '"«(?) cannot exceed its norm, hence it is bounded by (13.11 1. Thus, discarding a factor {2ny^^^ < 1, 
we have for all k and n 

s((l+a)eW)» 



\u(k, 01 < e-'"^l*l 
where e{K) is defined in ( 12.7b . Now choosing 



exp(ln J—e"'^"' - no-\k\j. 



Ut 



(3.6) 
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we obtain with n = In \k\ll3it,6 the following estimate 



\u{k, t)\ < exp 



- — \k\ In Ml — ) 

/iK,6 ^ o-ln\k\ I 



<e-^'*""'*'. (3.8) 



Remark 3.1 Since s and 6 can be chosen arbitrarily small and k arbitrarily large. Corollary 3.1 implies that, in 
terms of the dimensionless wavenumber k - 2crk, the Fourier amplitude has a bound (at high enough k) of the form 
^-ckink ^jjy (-^ ^ 1/(2 ln2). We shall see that the upper bound for the constant C can probably be improved to 
l/ln2. 

Proof of Theorem 3.1 We proceed again by induction. We assume that the following inequalities hold 



1/2 



|e-^-„(0f < |l, (3.9) 



and 



(s)fds<^, (3.10) 



where /r„ andfl„ > 1 are still to be determined. Starting from expNSE ( 12.3b . we take the inner product with e^'^"^'^°^'^'^^M. 
Then we obtain from Proposition 2. 1 with a - no- and p = 2cr 



2 dt 

< |e"™ "B(M,M)| |e' 



«-^'"b(«,«)| |e<'-2).A.'^ I < CA(4(y8)f ^' |e«-^"'«f"*'' |e(«+2).A./^ (3.11) 



< ^^-dWcf )(/.(;6))''^*^ (1 +a(.))^« |e-^"^«f''^' + ^ |e(-2)-^"^«f . 



4/f 

Then it follows that 

^ £^^-^«Cr (H^f le"-'^/^*" 

By using the induction assumption we obtain 



(3.12) 



||e(-i)-^"'«| +^|e<"^2KA"^^| <C"(^) , (3.13) 



where we have set 



C".^^-«Cr(/.(^)f*^(^| . (3.14) 
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Renaming the time variable in (13.13b from t to s' and integrating over s' from i to f (with < s < t < T) 



|2 



|e(''+')'^^"'«(f)f + // J' |e("+2K^"'M(.')| ds' 
< 1 C" Kf ^' ( \—\ + |e(«+i)-^"'M(.)|' 

< C'Kf ^'—^ + |e('-'»-^'%(*)|' . 

a„e{K) - 1 " I I 

Omitting the positive integral term on the left-hand side of the inequality we obtain 

|e("+')-^"'M(f)r < ^ C"Kf'^—^ + |e("+')-^"'«(^)f . 

I I fl„e(/f) - 1 " ,5«»eW-i I ^ '\ 

Choosing 1 <y < - {LjUT)^ and integrating over s from tjy to t we obtain 

yl ^'^l ane{K)-\a„e(K)-2 " \^ t > ^t' I 

J,/yl I fl„eW-lfl„eW-2 " W'' "^t' ' 

where we have used the induction assumption (13.101 1. We obtain thus the following estimate 

u(t)\ < C" Kf'H ^ + kJ ^ , 



^(n+l)o-A"-, 



which holds for every < f < T. 

To estimate e^"^^'°"'^ "m(.s) ds we integrate ( 13.13b from t to T: 



|2 



< C" ^ C" Kf '^—^ + |e<"+"'^^"'M(f)r . 

fl„eW-l " f«"ew-i 1 ''I 

Omitting the first term on the right-hand side and using (13.18b we obtain 

f |e(«+2KA"^«(,)|' ds < C" \ C" Kf '^^rn 
J, I I a„e(K) - 1 " f«»ew-i 

y- \ a„e{K)-lane{K)-2 " ^t' y - 1 t ' ' 

We conclude that since a„ > 1 and 

fl„ + 1 < a„{2 + 2j_^ ) ^ '^"^('*')' 

for a suitable constant £ > we have 

\a„e(K) 



14 
and 



Since, for t < T, 



it follows that 
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\a„e(K) 



~t~~r~ vufi 7 " vut) T 



and 



J|e«-'»»(.,|.£*:r(^) (I) 



This finishes the induction step. 

From the above follows that we can take 



- a„(l + 6) 



Ak-5 
2k- 5' 



(3.23) 



(3.24) 



(3.25) 



(3.26) 



Note that in the induction step we use the assumption that a,, > 1 . This fixes the value of uq - \. The solution of the 
recursion relations is given by 



fl„ = (1 +5) 



Ak-5 
2k- 5 



lnK„ = InE 



a„ - 1 



+ a„ InKo. 



(3.27) 



Finally, choosing a sufficiently large number A we get the desired estimates ( 13.11 ) and i3.2i . This concludes the proof 
of Theorem 3.1. 



4. Remarks and extensions for the main results 

Although our main theorems are stated for the case of the 3D expNSE, their statements and proofs are easily extended 
mutatis mutandis to arbitrary space dimensions d: with exponential dissipation for any d the solution is entire in the 
space variables and the decay of Fourier coefficients is bounded by exp(-C|^| In 1^1) for any C < C* = 1/(2 ln2). Some 
of the intermediate steps in the proof, such as the formulation of Agmon's inequality, change with d but not the result 
about the constant 1 / (2 In 2). 
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We can also easily change the functional form of the dissipationO One instance is a dissipation operator T) with a 
Fourier symbol e^""'*^'" with < a < 1 . One can prove that the solution in this case satisfies 

J]e2«W"|fi(fc,0l'<^ + C„, (4.1) 

for all t e (0, T] and for all n. Hence the solution in this case belongs to but is not necessarily an entire function. In 
fact it belongs to the Gevrey class Gija- Gevrey regularity with < a < 1 does not even imply analyticitylf] Actually, 
with such a dissipation, the solutions are analytic even when a < \. We shall return to this case of dual Gevrey 
regularity and analyticity in Section lSTTI 

Next, consider the case a > \. The dissipation has a lower bound of the ordinary exponential type, so that the entire 
character of the solution is easily established. However, for a > 1 the bound exp(-C|A;| In \k\) can be improved in its 
functional form, as we shall see in Section lSTI 

Obviously, the results of Sections |2] and [3] do still hold if we change the functional form of the Fourier symbol of 
the dissipation at low wavenumbers |^| while keeping its exponential growth at high wavenumber. One particularly 
interesting instance, to which we shall come back in the next Section on the Burgers equation and in the Conclusion, 
is "cosh dissipation", namely a Fourier symbol - cosh(A;/A;d)) with yu > 0. The dissipation rate at wavenumber 
much smaller than is then v\k\^ with v — yL(/(2A:d^), just as for the ordinary Navier-Stokes equation. 

It is worth mentioning that the key results of Sections |2] and [3] still hold when the problem is formulated in the 
whole space W' rather than with periodicity conditions. Similarly they should hold on the sphere 5^, a case for which 
spherical harmonics can be used (see 1*231). 

Of course the result on the entire character of the solution, when exponential dissipation is assumed, holds for a 
large class of partial differential equations. Besides the exponential modification of the Navier-Stokes equations it 
applies to similar modifications, for example, of the magnetohydrodynamical equations and of the complex Ginzburg- 
Landau equation 

du d^u T 

— -a—+l3u + y\u\\^Q, (4.2) 



Note that the proof of Proposition 2.1 and Theorem 2.1 holds mutatis mutandis if we replace, in the argument of the exponential, \k\ by a 

subadditive function of \k\ subject to some mild conditions, such as \k\'' with < a < 1. 

^ The special class when a = 1 of analytic functions is considered by some authors as one of the Gevrey classes 120i21] . 
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where 

Re a > , Re r > 0. (4.3) 

The main idea would be in proving the analogue of Proposition 2.1 for the corresponding nonlinear terms in the 
underlying equations following our proof combined with ideas presented in ||2T1| and ll24l . 



5. The case of the ID Burgers equation 

The (unforced) one-dimensional Burgers equation with modified dissipation reads: 

du du ^ _ 
— +u—^-luDu, (5.1) 

u(x, 0) = uo(x). (5.2) 

We shall mostly consider the case of the cosh Burgers equation when D has the Fourier symbol -/i(l - cosh(k/ki)). 
Since the cosh Burgers equation is much simpler than expNSE we can expect to obtain stronger results or, at least, 
good evidence in favor of stronger conjectures. 

Let us observe that the cosh Burgers equation can be rewritten in the complexified space of analytic functions of 
z ■- X + iy as 

+ u(z, f)^^^ = [u(z + i /h, t) + u(z-i /h, t) - 2u(z, t)] . (5.3) 
at ox 1 

This is the ordinary Burgers equation with the dissipative Laplacian replaced by its centered second-order finite dif- 
ference approximation, differences being taken in the pure imaginary direction with a mesh Ijki. 

As akeady stated. Corollary 2.1 on the entire character of the solution and Corollary 3.1 on the bound of the 
modulus of the Fourier coefficients by exp(-C|^| In |^|) for any C < 1 /(2 In 2) hold in the same form as for the expNSE. 
Of course, if the finite differences were taken in the real rather than in the pure imaginary direction, the solution would 
not be entire. Actually, ( 15.31 ) relates the values of the velocity on lines parallel to the real axis shifted by +l/^d in 
the imaginary direction. It thereby provides a kind of Jacob 's Ladder allowing us to climb to complex infinity in the 
imaginary direction. This can be used to show, at least heuristically, that the complexified velocity grows with the 
imaginary coordinate y as exp 
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Such a heuristic derivation turns out to be equivalent to another derivation by dominant balance which can be done 
on the Fourier-transformed equation, the latter being not limited to cosh dissipation. Section lSTl is devoted to Fourier 
space heuristics for different forms of the dissipation. For exponential and cosh dissipation this suggests a leading-order 
behavior of the Fourier coefficients for large wavenumber of the form exp(-C*|A:| In with C* - 1 / In 2, a substantial 
improvement over the rigorous bound. Various numerical and semi-numerical results, discussed in Section l572l and l531 
support this improved result. 



5.7. Heuristics: a dominant balance approach. We want to handle dissipation operators D with an arbitrary positive 
Fourier symbol, taken here to be e*^'*^ where G(k) is a real even function of the wavenumber k €Z which is increasing 
without bound for A: > 0. It is then best to rewrite the Burgers equation in terms of the Fourier coefficients. We set 

M(j(;,f) = ^e'*-'M(A:,0, (5.4) 



teZ 



and obtain from ( 15. It 



dii{k, t) i k 



p+cj=k 

This is the place where we begin our heuristic analysis of the large-wavenumber asymptotics. First, we drop the time 
derivative term since it will turn out not to be relevant to leading order (A suitable Galilean change of frame may be 
needed before this becomes true.) For simplicity we now drop the time variable completely. The next heuristic step is 
to balance the moduli of the two remaining terms, taking 

\um ~ e"''^*', (5.6) 

where F{k) is still to be determined but assumed sufficiently smooth and the symbol ~ is used here to connect two 
functions "asymptotically equal up to constants and algebraic prefactors" (in other words, asymptotic equality of 
the logarithms). The convolution in ( 15.51 ) can be approximated for large wavenumbers by a continuous wavenumber 
integral ~ / t-'^^P^-^^'^-P^dp. Next we evaluate the integral by steepest descent, assuming that the leading order comes 
from the critical point p - k/2, where the /j-derivative of F(p) + F{k - p) obviously vanishes. This will require that 
this point be truly a minimum of F{p) + F{k - p). Balancing the logarithms of the nonUnear term and of the dissipative 
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term we obtain the following simple equation for the function F(ky. 

2F^'^'j = F(k)-G(k). (5.7) 

This is a linear first order finite difference equation (in the variable In A:) which is easily solved for values of the 
wavenumber of the form k -2": 



F(2") = 2" 



2" 



(5.8) 



For exponential dissipation (and for cosh dissipation when \k\lki » 1), we have G{k) ^ 2o-\k\ and we obtain from 
( 15.81 1. to leading order for large positive k 

F{k) ^ -^k\nk; k:=2crk=^. (5.9) 
In 2 % 



If this heuristic result is correct - and the supporting numerical evidence is strong as we shall see in Sections [5.21 and 
15.31 - the estimate given by Corollary 3.1 (adapted to the Burgers case) that \u{k)\ < Q-^Wl'^niki ^^j. sufficiently large |^| 
and any C < C* = 1 /(2 In 2) still leaves room for improvement as to the value of C* . It can be shown that this dominant 
balance argument remains unchanged if we reinsert the time-derivative term, since its contribution is easily checked 
to be subdominant. Actually, the conjecture that the solution of NSE is entire with exponential or cosh dissipation 
was based on precisely this kind of dominant balance argument, which suggests a faster-than-exponential decay of the 
Fourier coefficients. 

When G{k) - 2o-\k\" with o' > 1 we obtain to leading order 



.0 



This is an even faster decay of the Fourier coefficients than in the exponential case (I2.3l ll' 

It is easily checked that for a > I the condition of having a minimum of F(p) + F(k - p) at p - k/2 is satisfied. 
If however we were to use ( |5.10t for < a < 1 the condition would not be satisfied. In this case it is easily shown 
for the Burgers equation and the NSE, by using a variant of the theory presented in Section |2] that the solution is in 
the Gevrey class Gi/q. in the whole space R"'. It is actually not difficult to show that the solution is also analytic when 



Actually, one can show, for the Burgers equation and the NSE that when a > 1 the Fourier coefficients of the solution decay faster than 
expi-Qkl"-^) for any e > 0. 
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< a < 1, in a finite strip in C' about the real space R''. For this it suffices to adapt to the proof of analyticity given 
for the ordinary NSE under the condition of some mild regularity. Such regularity is trivially satisfied with the much 
stronger dissipation assumed here 1201 P We also found strong numerical evidence for analyticity. 

It is of interest to point out that, although analyticity is a stronger regularity than Gevrey when < a < 1, the 
Gevrey result implies a decay of the form exp(-C|A;|''ln independently of the viscosity coefficient jj., whereas 
analyticity in a finite strip gives a decay of the form exp(-77|^|) where 77 depends on fi ||26l . 



5.2. Spectral simulation for the Burgers case. Here we begin our numerical tests on the ID Burgers equation. So 
far we have a significant gap in the value of the constant C appearing in the e"'''*^' estimate of Fourier coefficient, 
between the bounds and a heuristic derivation of the asymptotic behavior In this section we shall exclusively consider 
the case of the unforced Burgers equation with initial condition U()(x) - - sin x and dissipation with a rate 1 - cosh k. 
(Thus, n - I and - 1.) The numerical method is however very easily extended to other functional forms of the 
dissipation and other initial conditions. The spectral method is actually quite versatile. Its main drawback will be 
discussed at the end of this section. 

The standard way of obtaining a high-orders scheme when numerically integrating PDE's with (spatial) periodic 
boundary conditions is by the (pseudo)-spectral technique with the 2/3 rule of alias removal ||27ll . The usual reason 
this is more precise than finite differences is that the truncation errors resulting from the use of a finite number of 
collocation points (and thus a finite number of Fourier modes) decreases exponentially with if the solution is 
analytic in a strip of width 6 around the real axis. Indeed this implies a bound for the Fourier coefficients at high |^| of 
the form \u{k)\ < e"*-'*' for any C < 5. In the present case, the solution being entire, the bound is even better. 

There are of course sources of error other than spatial Fourier truncation, namely rounding errors and temporal 
discretization errors. Temporal discretization is a non-trivial problem here because the dissipation grows exponentially 
with \k\ and thus the characteristic time scale of high-|^| modes can become exceedingly small. Fortunately, these 
modes are basically slaved to the input stemming from nonlinear interaction of lower-lying modes. It is possible to 
take advantage of this to use a slaving technique which bypasses the stiffness of the equation (a simple instance of this 

' The first results on analyticity, derived in the more complex setting of flow with boundaries, were obtained in 1251 . 
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Fig. 1. Discrepancy Discr, as predicted by 15.1 U . vs wavenumber k for data obtained by spectral simulation in double precision; rounding noise 
becomes significant beyond wavenumber 17. 



phenomenon is described in Appendix B of f28l). We use here the slaved scheme Exponential Time Difference Runge 
Kutta 4 (ETDRK4) of EH with a time step of 10"^ Q 

As to the rounding noise, it is essential to use at least double precision since otherwise the faster-than-exponential 
decrease of the Fourier coefficients would be swamped by rounding noise beyond a rather modest wavenumber Even 
with double precision, rounding noise problems start around wavenumber 17, as we shall see. Hence it makes no sense 
to use more than, say, 64 collocation points, as we have done. 

Fig.[T]shows the discrepancy 

ln|M(;t, 1)1 1 



Discr (/t) 



(5.11) 



|A:|ln|;t| hi 2 

which, according to heuristic asymptotic theory ( |5.6I )-( |S!9] |, should converge to as \k\ oo. 

It is seen that the discrepancy falls to about 3.5% of the nominal value 1/ In 2 before getting swamped by rounding 
noise around wavenumber 17. 



* This is far larger than would have have been permitted without the slaving. Actually it can still be increased somewhat to 5 X 10 ^ without 
affecting the results. 
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k 

Fig. 2. Discrepancy Discr vs wavenumber k for the same data as in Fig.[T] but processed by a 5-stage asymptotic extrapolation method. 

It is actually possible to significantly decrease the discrepancy by using a better processing of the numerical output, 
called asymptotic extrapolation, developed recently by van der Hoeven [30] and which is related to the theory of 
transseries II31II32I . The basic idea is to perform on the data a sequence of transformations which successively strip 
off the leading and subleading terms in the asymptotic expansion (here for large \k\). Eventually, the transformed data 
allow a very simple interpolation (mostly by a constant). The procedure can be carried out until the transformed data 
become swamped by rounding noise or display lack of asymptoticity, whichever occurs first. After the interpolation 
stage, the successive transformations are undone. This determines the asymptotic expansion of the data up to a certain 
order of subdominant terms. An elementary introduction to this method may be found in ||33l . from which we shall 
also borrow the notation for the various transformations: I for "inverse", R for "ratio", SR for "second ratio", D for 
"difference" and Log for "logarithm". The choice of the successive transformations is dictated by various tests which 
roughly allow to find into which broad asymptotic class the data and their transformed versions fall. In the present case, 
the appropriate sequence of transformations is: Log, D, D, I, D. Because of the relatively low precision of the data it 
is not possible to perform more than five transformations, so that the method gives us access only to the leading-order 
asymptotic behavior, namely \u{k, 1)| ^ e-c,l*|in|<:| jj j^^^y shown that the constant C* - where w^^' is the 
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constant value of the high-|A;| interpolation u^^\\k\) after the 5th stage of transformation. Fig.|2]shows the discrepancy 
u^^\\k\) + In 2 . The absolute value of the discrepancy lingers around 0.002 to 0.005 before being swamped by rounding 
noise at wavenumber 17. Thus with asymptotic extrapolation the discrepancy does not exceed 0.7% of the nominal 
value In 2. The accuracy of the determination has thus improved by about a factor 5, compared to the naive method 
without asymptotic extrapolation. 

To improve further on this result and get some indication as to the type of subdominant corrections present in 
the large-wavenumber expansion of the Fourier coefficients it would not suffice to increase the spatial resolution, 
since rounding noise would still swamp the signal beyond a wavenumber of roughly 17. Higher precision spectral 
calculations are doable but not very simple because high-precision fast Fourier transform packages are still in the 
experimental phase. 

In the next Section we shall present an alternative method, significantly less versatile as to the choice of the initial 
condition because it exploits the algebraic structure of a certain special class of solutions, but which also allows to 
work easily in arbitrary precision and thus to make better use of asymptotic extrapolation for determining the constant 
C*. 

5.3. Half-space (Fourier) supported initial conditions. So far we have limited ourselves to initial conditions that are 
real entire functions. Hence the Fourier coefficients had Hermitian symmetry: uoi-k) - Mg(A:), where the star denotes 
complex conjugation. With complex initial data there are no analytical results when the dissipation is exponential, even 
when the initial conditions are entire because the energy conservation relation - now about a complex-valued quantity 
- ceases to give L^-type bounds. Actually, as already pointed out, Li and Sinai ifTSl showed that the 3D NSE can display 
finite-time blow up with suitable complex initial data. It is however straightforward to adapt to complex solutions the 
heuristic argument of Section ISTT] and to predict a high-wavenumber leading-order term exactly of the same form as 
for real solutions. This is of interest since we shall see that there is a class of periodic complex initial conditions for 
which, provided the Burgers equation is written in terms as the Fourier coefficients as in 1341 and liT3l , any given 
Fourier coefficient can be calculated at arbitrary times f with a finite number of operations, most easily performed on 
a computer, by using either symbolic manipulations or arbitrary-high precision floating point calculations. 
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For the case of the Burgers equation, this class consists of initial conditions having the Fourier coefficients sup- 
ported in the half line k > OO We shall refer to this class of initial data as "half-space (Fourier) supported". 

Because of the convolution structure of the nonlinearity when written in terms of Fourier coefficients, it is obvious 
that with an initial condition supported in the k > half line, the solution will also be supported in this half line. A 
similar idea has been used in three dimensions for studying the singularities for complex solutions of the 3D Euler 
equations ll35l . 

Specifically, we consider again the ID Burgers equation (I5.1l l-( |S!2] | with 27r-periodic boundary conditions for t > 
0, rewritten as ( 15. St , in terms of the Fourier coefficients u{k, f), assumed here to exist. The ^-dependent real, even, 
non-negative dissipation coefficient yue*^*** is denoted by p(k). The initial conditions uoik) (k = 1,2, . . .) are chosen 
arbitrarily, real or complex. We then have the following Proposition, which is of purely algebraic nature: 

Proposition 5.1 Eq. i5.5\ with the initial conditions u(k, 0) — (k < 0) and u{k, 0) = uoik) (k — 1,2,...) defines, 
for all k > and t > 0, u(k, t) as a polynomial functions of the set ofuo(p) with < p < k. 

Proof From ( 15.51 1. after integration of the dissipative term, we obtain, for A: > and t > 

■ , k-l 

u(k,t)^e-'P"'^uo(k)-— I dse-^'-'^P^'^'^y u{p,s)u(k- p,s). (5.12) 

Observe that ii(k, t), given by ( 15.12b . involves u(k, 0) (linearly) and the set of Fourier coefficients u(p, s) for I < p < 
k-l and < s < t (quadratically). The proof follows by recursive use of this property for k, k- 1, . . . , iQ 

Note that the solution can be obtained without any truncation error on a computer, using symbolic manipulation. 
Alternatively it can be calculated in arbitrary high-precision floating point arithmetic. 

Now we specialize the initial condition even further, by assuming that the only Fourier harmonic present in the 
initial condition has ^ = 1 111 We take 



uo(x)^iJie\ (5.13) 
for which mo(1) = i^t, while all the other coefficients vanish. 



' If the coefficient for wavenumber k = Ois non- vanishing a simple Galilean transformation can be used to make it vanish. 

This proposition has an obvious counterpart for the NSE in any dimension when the the Fourier coefficients of the initial condition are compactly 

supported in a product of half-spaces. 
' ' What follows can be easily extended to the case of a finite number of non- vanishing initial Fourier harmonics. 
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Fig. 3. Fifth stage of asymptotic inteipolation showing the discrepancy with respect to the heuristic prediction C* = 1 / In 2. Note the decaying 
oscillatory behavior. 



Setting, for ^ > 0, 



v(k, f) := 



iiik, f) 



(5.14) 



we obtain from ( I5.12l i by working out the power series to the second order the following fully expUcit expressions of 
the first two Fourier coefficient at any time t >0: 



v(l,f) = e 



-p{\)t 



y(2,f) 



-2p(l)f 



p(2)-2p(l) p(2)-2p(l) 



(5.15) 



For the one-mode initial condition with and p{k) - e'*^ (exponential dissipation with - \ and {k^ - 1), we 

have calculated the Fourier coefficients vik, 1) for k - 1, 2, . . . , 24 using Maple symbolic calculation with a forty-digit 
accuracy. 

The data have then been processed using asymptotic extrapolation as in Section l572l with the same five transfor- 
mations Log, D, D, I, D. Fig. |3]shows the discrepancy v^^\k, 1) -i- In 2 between the 5th stage of interpolation and the 
prediction from the dominant balance argument. It is seen that the discrepancy drops to -2.3 x 10"^. Thus the relative 
error is about 3 x 10"^. The oscillations in the discrepancy, if they continue to higher wavenumbers would indicate 
that the first subdominant correction to the asymptotic behavior of the Fourier coefficient e"('/'"2)|*in|i| ^ prefactor 
involving a complex power of the wavenumber. 
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Finally, we address the issue of what kind of solution we have constructed by this Fourier-based algebraic method. 
Is it a "classical" sufficiently smooth global-in-time solution of the Burgers equations ( |5.U -( IJ!2] i written in the physical 
space? We have here obtained strong numerical evidence that the Fourier coefficients decrease faster than exponentially 
with the wavenumber and thus define a classical solution which is an entire function of the space variable. This is 
however just a conjecture. The tools used in Section|2]to prove the entire character of the solution rely heavily on the 

u 

definite positive character of the energy, a property lost with complex solutions! n 



6. Conclusions 

In this paper we have proved that for a large class of evolution PDE's, including the 3D NSE, exponential or faster- 
growing dissipation implies that the solution becomes and remains an entire function in the space variables at all 
times. Exponential growth constitutes a threshold: subexponential growth with a Fourier symbol e'*^'", where < a < 1 
makes the solution analytic (but not entire) as is the case in 2D (and generally conjectured in 3D). Furthermore, for 
the 3D NSE and the ID Burgers equation with a dissipation having the Fourier symbol //e'*^'^*^'', we have shown that 
the amplitude of the Fourier coefficients is bounded by e^c\k\in\k\ (y^jjgj-g ^ k/k^) for any C < 1/(2 In 2). For the 
case of the ID Burgers equation we have good evidence that this can be improved to C < C* = l/ln2 since the 
high-|^| asymptotics seems to have a leading term precisely of the form Q-c*\k\\n\k\.^ evidence comes both from a 
heuristic dominant balance argument and from high-precision simulations. The heuristic argument can actually be 
carried over somewhat loosely to the expNSE in any dimension: again the dominant nonlinear interaction contributing 
to wave vector k comes from the wave vectors p - q - k/2; actually the condition of incompressibility kills nonlinear 
interactions between exactly parallel wave vectors but this is only expected to modify algebraic prefactors in front of 
the exponential term. 

We thus conjecture that the C* = l/ln2 also holds for expNSE in any space dimension c/ > 2. Of course there 
is a substantial gap between the bound and the conjectured asymptotic behavior It seems that such a gap is hard to 
avoid when using L^-type norms. For proving the entire character of the solution such norms were appropriate. Beyond 
this, it is appears more advisable to try bounding directly the moduli of Fourier coefficients by using the power series 

It is however not difficult to prove, for short times, that our solution is also a classical entire solution. 
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method 013II34I . A first step in this direction would be to prove that C* = 1/ In 2 for initial conditions whose Fourier 
coefficients are compactly supported in a product of half-spaces, of the kind considered in Section [531 '-^ 

Exponential dissipation differs from ordinary dissipation (with a Laplacian or a power thereof) not only by giving 
a faster decay of the Fourier coefficients but by doing so in a universal way: with ordinary dissipation the decay of the 
Fourier coefficient is generally conjectured to be, to leading order, of the form e where 77 depends on the viscosity v 
and on the energy input or on the size of the initial velocity; with exponential dissipation the decay is Q-c,\k\in\k\ ^vhere 
C* - 1 / In 2 and thus depends neither on the coefficient which plays the role of the viscosity nor on the initial dataB 
As a consequence, it is expected that exponential dissipation will not exhibit the phenomenon of dissipation-range 
intermittency, which for the usual dissipation can be traced back either to the fluctuations of 77 1 36 1 or to complex 
singularities of a velocity field that is analytic but not entire ll37l . 

Finally some comments on the practical relevance of modified dissipation. First, let us comment on "hyperviscos- 
ity", the replacement of the (negative) Laplacian by its power of order a > 1 . Of course we know that specialists of 
PDE's have traditionally been interested in the hyperviscous 3D NSE, perhaps to overcome the frustration of not being 
able to prove much about the ordinary 3D NSE. But scientists doing numerical simulations of the NSE, say, for engi- 
neering, astrophysical or geophysical applications, have also been using hyperviscosity because it is often believed to 
allow effectively higher Reynolds numbers without the need to increase spatial resolution. Recently, three of us (UF, 
WP, SSR) and other coauthors have shown that when using a high power a of the Laplacian in the dissipative term for 
3D NSE or ID Burgers, one comes very close to a Galerkin truncation of Euler or inviscid Burgers, respectively |38|. 
This produces a range of nearly thermalized modes which shows up in large-Reynolds number spectral simulations as 
a huge bottleneck in the Fourier amplitudes between the inertial range and the far dissipation range. Since the bottle- 
neck generates a fairly large eddy viscosity, the hyperviscosity procedure with large a actually decreases the effective 
Reynolds number. 

Next, consider exponential dissipation. In 1996 Achim Wirth noticed that when used in the ID Burgers equation, 
cosh dissipation produces almost no bottleneck although it grows much faster than a power of the wavenumber at high 

Progress on this issue has been made and will be reported elsewhere. 

It does however depend on the type of nonlinearity. For example with a cubic nonlinearity the same kind of heuristics as presented in Section lsTT] 
predicts a constant C* = 1 / In 3. 
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wavenumbers 1391 . It is now clear that such a dissipation will produce a faster-than-exponential decay at the highest 
wavenumbers. But at wavenumbers such that |^| <k A:d a dissipation rate -jt/(l - coshk/k^) reduces to jj\k\^ /(2kd'^), to 
leading order, which is the ordinary (Laplacian) dissipation. With the ordinary ID Burgers equation it may be shown 
analytically that there is no bottleneck. For the ordinary 3D NSE, experimental and numerical results show the pres- 
ence of a rather modest bottleneck (for example the "compensated" three-dimensional energy spectrum |A:|^''^''£'(|^|) 
overshoots by about 20%.). If in a simulation with cosh dissipation // and are adjusted in such a way that dissipation 
starts acting at wavenumbers slighter smaller than k^, the beginning of the dissipation range will be mostly as with an 
ordinary Laplacian, that is with no or little bottleneck!^ At higher wavenumbers, where the exponential growth of the 
dissipation rate is felt, faster than exponential decay will be observed. In principle this can be used to avoid wasting 
resolution without developing a serious bottleneck. Faster than exponentially growing dissipation, e.g. jU ^e*^!*^'^*'''^ ~ 
may be even better because the prediction is that the Fourier coefficients will display Gaussian decayl^ Testing the ad- 
vantages and drawbacks of different types of faster-than-algebraically growing dissipations for numerical simulations 
is left for future work. 
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